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ABSTRACT 


Carotenoids, which generate yellow, orange, and red 
colors, are crucial pigments in avian plumage. 
Investigations into genes associated with carotenoid- 
based coloration in avian species are important; 
however, such research is difficult because 
carotenoids cannot be synthetized in vertebrates as 
they are only derived from dietary sources. Here, the 
golden pheasant (Chrysolophus pictus) was used as 
a model in analysis of candidate gene expression 
profiles implicated in carotenoid binding and 
deposition. Using mass and Raman spectrometry to 
confirm the presence of carotenoids in golden 
pheasant feathers, we found C4oH540 and C4oHs5602 
in feathers with yellow to red colors, and in the 
rachis of iridescent feathers. The global gene 
expression profiles in golden pheasant skins were 
analyzed by RNA-seq and all six carotenoid binding 
candidate genes sequenced were studied by real- 
time PCR. StAR4, GSTA2, Scarb1, and APOD in 
feather follicles showed different expressions in red 
breast and orange nape feathers compared with that 
of iridescent mantle feathers. Further comparison of 
golden pheasant yellow rump and Lady Amherst's 
pheasant (Chrysolophus amherstiae) white nape 
feathers suggested that GSTA2 and APOD played a 
potential role in carotenoid-based coloration in 
golden pheasant. 
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INTRODUCTION 


Golden pheasant (Chrysolophus pictus) is a world-renowned 
avian species with particularly attractive and varied plumage. As 
one of the most prevalent integument pigments in birds, 
carotenoids are usually deposited in bare body parts and 
produce yellow, orange, and red coloration in feathers (usually 
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in the distal parts). To date, around 30 different carotenoids 
have been identified in hundreds of avian species (Hill & 
McGraw, 2006; LaFountain et al., 2010). Furthermore, this class 
of pigment has been confirmed in 95 of 236 extant bird families, 
including three species of Galliformes, such as the golden 
pheasant (Thomas et al., 2014a). Unlike other pigments, such 
as melanin, carotenoids are difficult to study in regards to the 
genetics of bird coloration because knowledge of the molecular 
genetics of carotenoid-based coloration is incomplete (Roulin & 
Ducrest, 2013). Under this limitation, transportation, deposition, 
and metabolism-related proteins and their encoding genes have 
been studied in consideration of the dietary-derived characteristics 
of these pigments (Pointer et al., 2012; Walsh et al., 2012). In 
chicken skins, mutations of beta-carotene dioxygenase 2 
(BCDO2) inhibit gene expression and result in carotenoid 
accumulation (Eriksson et al., 2008). In bird feathers, 
carotenoids are bound directly to keratin or other proteins 
(Stradi et al., 1995; Mendes-Pinto et al., 2012), so it is generally 
accepted that carotenoid-binding protein (CBP) is the starting 
point to investigate carotenoid-based coloration of plumage at 
the molecular level. 

A previous report reviewed 11 plausible candidate genes 
associated with carotenoid uptake, binding, and deposition, 
including steroidogenic acute regulatory protein 4 (SfAR4), 
glutathione S-transferase alpha 2 (GSTA2), StAR-related lipid 
transfer domain containing 3 (STARD3, also known as MLN64), 
scavenger receptor class B, member 1 (Scarb1, also known as 
SR-B/), apolipoprotein D (APOD), and lipid storage droplet 2- 
like (PLIN) genes (Walsh et al., 2012). As a member of the 
StAR protein family, StAR4 specifically binds lutein with high 
affinity, and is a possible candidate for CBP in human retina 
(Bhosale et al., 2009). Scarb1 is involved in carotenoid 
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transport into retinal pigment epithelium (RPE) cells, and gene 
silencing can result in a marked inhibition of carotenoid uptake 
(During et al., 2008). Although APOD is a carotenoid-binding 
gene in crustaceans and mammals, a paralogous gene has 
since been characterized in chickens, and is expressed in skins 
and developing feathers (Ganfornina et al., 2005). GSTAZ2 is the 
closest avian homologue of the carotenoid binding gene 
GSTP1 in mammalian retina and is involved in plumage 
dichromatism (Walsh et al., 2012; Zhang et al., 2014). It is 
predicted, therefore, that these genes might be highly 
expressed in carotenoid feathers. 

Prior to pigment identification, we analyzed the expression 
levels of six carotenoid deposition genes in golden pheasant 
feather follicles, in which the corresponding feathers may or 
may not have contained carotenoids. 


MATERIALS AND METHODS 


Animals and samples 

Three male artificially bred golden pheasants were used in this 
study. The birds were provided by the Yuanfeng Wildlife Farm 
(Jilin, China) and were housed in cages with free access to 
water and commercial food. This study was approved by the 
Institutional Animal Care and Use Committee of the Inner 
Mongolia University. For RNA extraction, skins and feather 
follicles were sampled. Briefly, feathers of various colors were 
plucked from the breast, nape, rump, and mantle of the body 
(Figure 1). Regrown feather shafts and skins were then 
collected and stored in liquid nitrogen (Poelstra et al., 2014). 
For pigment extraction, feathers were plucked from 
corresponding regions directly, then placed in envelopes and 
stored in the dark. Blood samples were drawn from the neck 
veins using vacuum blood collection tubes, with the samples 
then centrifuged in heparinized microcapillary tubes for 10 min 
at 1 500 g. The plasma was then removed and stored at —80 °C 
in 1.5 mL Eppendorf tubes until further analysis (McGraw et al., 
2002a). 


Nape 
Mantle 





Breast 


Figure 1 Feathers and feather follicles sampled from the breast, 
nape, mantle and rump of the golden pheasant 

Feathers from breast, nape and rump were red, orange and yellow in 
the distal parts, while feathers from the mantle were iridescent green. 
Feather follicles were sampled during the growing period. 


RNA-seq 
The isolated mRNA from golden pheasant skins was 


fragmented and used for cDNA synthesis. The short cDNA 
fragments were treated following the manufacturer’s 
instructions and sequenced using an Illumina Hiseq 2000 
(Illumina, San Diego, USA). The produced reads were 
assembled by mapping to the golden pheasant genome we 
sequenced and annotated based on homology prediction. The 
expression abundance of each assembled transcript was 
measured through reads per kilobase transcriptome per million 
mapped reads (RPKM) (Mortazavi et al., 2008). The 
expressions of each read between sample pairs were 
calculated according to RPKM values in different samples. A 
minimum two-fold difference was considered as an expression 
difference with a false discovery rate (FDR) $0.05 (Benjamini et 
al., 2001). GO enrichment analysis was performed using GO: 
TermFinder software (http: //smd.stanford.edu/help/ GO- 
TermFinder/GO_TermFinder_help.shtml/). 


Carotenoids extraction 

Feather and plasma carotenoid pigments were extracted, 
respectively, using previously published methods (McGraw 
et al., 2002b, 2005). In brief, feathers from various parts of 
the body were soaked in ethanol and hexane, in turn, in 
order to remove surface lipids. Feather vanes were 
trimmed into pieces and placed into 10 mL glass tubes. 
After that, 1 mL of acidified pyridine was infused to cover 
the pieces. The tubes were then filled with argon gas, 
capped tightly, and placed in a 95 °C water bath for 3 h 
until the solution became colorful. After the thermochemical 
procedures, the cooled tubes were added with 2 mL of 
pure water and 1 mL of hexane: tert-butyl methyl ether (1:1) to 
separate the carotenoid pigments from the solution. The 
mixtures were shaken strongly for 2 min and centrifuged at 
3000 r/min for 5 min, with the supernatants then 
transferred to clean glass tubes, dried under nitrogen gas, 
and stored at —80 °C until use. For blood samples, the 
extraction procedures were the same as above (McGraw et 
al., 2002a). Ten microliters of plasma with 100 uL of 
ethanol were mixed in 1.5 mL Eppendorf tubes, followed by 
the addition of 100 uL of tert-butyl methyl ether. After being 
vortexed, the tubes were centrifuged at 3 000 r/min for 3 
min, with the same treatment as that for plumage pigment 
extraction. 


Mass spectrometry 

The extracted pigments were used for molecular weight 
analysis using an HPLC-MS instrument composed of a LC- 
20AD system (Shimadzu, Japan) fitted with a Waters™ 
C18 column (5 um, 150 mmx2.1 mm) and an LTQ Velos 
Orbitrap (Thermo, Germany). Samples were dissolved in 
methanol and injected into the instrument at 0.2 mL/min 
and analyzed in the positive mode with the capillary heated 
to 350 °C and at a cone voltage of 4.5 kV. A gradient 
solvent system was used for compound elution, with the 
percentage of acetonitrile varying from 5% to 95% 
corresponding to the methanoic acid solution (v=5%) in the 
first 11 min, before decreasing to 5% in the next 3 min, and 
finally being held isocratically for 1 min. 
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Raman spectroscopy 

Feather samples were examined in a Labram HR1800 
spectrometer (HORIBA Jobin Yvon, France) without pretreatment 
using an excitation laser source at 514 nm and 0.3 mW power. 
Single spectra were obtained using a 100x confocal objective 
and a grating of 600 lines/mm. Data were analyzed by LabSpec 5 
software. Standards of lutein (Chromadex, USA) and zeaxanthin 
(Sigma, USA) were used in comparative analysis. 


Real-time PCR 
Sequences of candidate genes were obtained from the RNA- 
seq data. The predicted sequences were used to design 


Table 1 Primers used for real-time PCR 


Gene Forward primer (5'-3') 
StAR4 CAAGCTATGAAGATGGGCTTC 
GSTA2 TCCCTTTTCAAGCAGCCGAT 
STARD3 GACCACGCACAGCCTGAA 
Scarb1 ACTTCTACAATGCTGACCCAA 
APOD GCGTCCGCTTCAACT GGT 
PLIN CCATCCAAAGTGCCAAGAG 
B-actin CTCCCTGATGGTCAAGTCAT 
RESULTS 


Carotenoid identification in golden pheasant plumage 
Pigments were separated from golden pheasant plasma and 
feathers according to solubility. Both procedures result in 
carotenoid transfer into the upper organic phase and profile 
colors in it if this kind of pigment is present. We found yellow 
colors emerged in the mentioned phases after yellow rump, 
orange nape, and red breast feathers were treated. This 
indicated that yellow fat-soluble pigments were contained in 
these feathers, but not in the iridescent mantle feathers. The 
same phenomenon occurred in the plasma samples. To confirm 
the molecular weight of the compounds, the extracted pigments 
from plasma and plumage were analyzed by mass 
spectrometry. Search results in databases (METACYC, LIPID, 
KEGG) showed that C4oHs40 and C4oHs602 were present in all 
samples. According to previous publication, lutein and 
zeaxanthin are isomers with the formula C4oHs5s602, and 
echinenone, which is a common carotenoid found in plumage, 
IS C4oHs40 (Hill & McGraw, 2006). 

A single iridescent feather taken from the mantle region of 
the golden pheasant was examined by Raman spectroscopy 
for carotenoid identification, as described previously 
(Thomas et al., 2014a, b). We took two collection points at 


Table 2 Statistics of transcriptome sequencing data 


primers on exon regions spanning introns using Primer Premier 
5 (Table 1). Housekeeping gene fB-actin was used for 
normalization. Total RNA was extracted from the feather 
follicle samples taken from each of the three birds using 
RNAiso Plus reagent (Takara, Japan). cDNA preparation with 
quantitative real-time PCR for each candidate gene was 
performed using a PrimeScript™ RT Reagent Kit with gDNA 
Eraser (Perfect Real Time) (Takara, Japan). Each sample was 
pooled with 5-10 feather follicles. Reactions were performed 
in an ABI PRISM 7500 Real-Time PCR System (Applied 
Biosystems, USA) with default parameters. Three replicates 
were used for each sample. 


Reverse primer (3'-5') TM (°C) 
ATCTGATTGCAAGGAATGCAC 60 
GCTGCCAGGCTGCAAGAAT 60 
AGGCTCTGGTGGATGAGGTA 60 
AGCTTGATGGAGCAATTCAT 60 
CGTGGGCATCATCTTGTCG 60 
CAGGTCTGCTTGGGCTTC 60 
TGGATACCACAGGACTCCAT 60 


the rachis and barb, respectively, and found characteristic 
bands of carotenoids in the rachis at 1 500-1 535 /cm 
(identified as v[C=C]), 1 145-1 165 /cm (identified as v[C-C}), 
and 1 000-1 010 /cm (identified as 6 [CH2]) (Thomas et al., 
2014b). For the barb, bands of eumelanin (Galvan et al., 
2013) were detected (Figure 2). 


Transcriptome sequencing and differential gene expression 
profiling 

Four libraries were constructed for golden pheasant skin 
transcriptome study, including the breast, nape, mantle, and 
rump. Raw data were sequenced by an Illumina Hiseq 2000 
and filtered to produce a total of 23.83 Gb clean reads. We 
mapped the reads to the golden pheasant genome to assemble 
the transcripts for each library. A summary of the mapping 
results is shown in Table 2. Transcriptome data were deposited 
in the Sequence Read Archive (SRA) database under 
accession no. SRP075618. 

A total of 50 228 expressed genes were predicted. To find 
differences in gene expression between two kinds of color 
morphs, we compared the RPKM values of the transcripts in 
mantle skin to the other three parts. Results showed that 452, 748, 
and 532 genes were upregulated and 686, 1 185, and 710 
genes were downregulated in the breast, nape and rump, 


Breast 
Total reads 68 796 446 
Total mapped reads 48 578 519 
Total mapped rate (%) 70.61 
Unique mapped reads 46 213 276 
Unique mapped rate (%) 67.17 
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Nape Mantle Rump 
64 410 352 64 910 570 66 653 538 
45 006 710 44 868 949 46 045 821 
69.87 69.12 69.08 
43 605 488 42 886 126 44 029 043 
67.70 66.07 66.06 
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Figure 2 Detected points (above) and corresponding spectrograms (below) by Raman spectroscopy in an iridescent mantle feather of 


the golden pheasant 


Purified standards of lutein and zeaxanthin were used as the positive controls. 


respectively. The genes with significantly differential expression 
were enriched in 46 pathways, including the retinol metabolism 
pathway, which might be involved in carotenoid metabolism 
(Table S1). We enriched and overlapped the genes that showed 
preferential expression in carotenoid-containing tissues (Figure 3). 
Sixty-six genes were significantly expressed in all three 
carotenoid-containing tissues, of which 57 were functionally 
annotated (Table S2). 


Expression of candidate genes in golden pheasant 
feather follicles 

To confirm the expressions of carotenoid binding related genes 
in feather follicles, six homologous candidate genes were 
determined using real-time PCR. Expression differences in 
developing feather follicles taken from the rump, breast, nape, 
and mantle regions of the male golden pheasant were 
determined and the expression patterns were compared with 


RNA-seq data (Figure 4). All candidate genes were expressed 
in all tissue samples. Comparison of the expression patterns of 
the six genes revealed no differences (P>0.05) between yellow 
rump and iridescent mantle feather follicles. Compared with 
mantle tissue, there was no difference (P>0.05) in the 
expression levels of STARD3 in the breast, nape and rump 
samples, and the expression of PLIN showed significant 
differences only in red breast feather follicles (P<0.05). The 
remaining five genes were differentially expressed in red and 
orange plumage. For red breast plumage, expressions of 
StAR4 and PLIN showed significant differences (P<0.05), while 
expressions of GSTA2, Scarb1, and APOD showed highly 
significant differences (P<0.01). For orange nape plumage, 
expressions of St(AR4 and APOD showed significant differences 
(P<0.05), while expressions of GSTA2 and Scarb1 showed 
highly significant differences (P<0.01). We also compared the 
expression levels of the candidate genes between golden 
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Figure 3 Analysis of differentially expressed genes (DEGs) upregulated in golden pheasant breast, nape, and rump compared with that 


in the mantle 


A-C: GO enrichment of DEGs in breast (A), nape (B), and rump (C). X axes indicate -Log2(P-value). D: Overlap of DEGs. 
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Figure 4 Expression levels of carotenoid candidate genes in different feather follicles of the golden pheasant (mean+tSD) by real- 


time PCR 


Columns indicate real-time PCR data, and lines indicate RNA-seg data. RPKM: Reads per kilobase transcriptome per million mapped reads. 
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pheasant yellow rump feathers and Lady Amherst's pheasantnape 
feathers, which are white with black stripes without carotenoid- 
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based coloration. Data showed the expressions of GSTA2 and 
APOD were significantly different (Figure 5). 
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Figure 5 Expression levels using real-time PCR of carotenoid candidate genes in different feather follicles between golden pheasant (G) 


rump and Lady Amherst's pheasant (L) nape 


DISCUSSION 


Thermochemical techniques for carotenoid extraction are based 
on weakening the non-covalent hydrogen bonds that bind 
pigments to proteins, thereby subsequently releasing pigments 
into solution (McGraw et al., 2005). In our study on carotenoid 
separation from golden pheasant plumage, we found that the 
organic phases of red, orange, and yellow tested barbs were 
pigmented after treatment. This phenomenon suggested the 
presence of lipochrome in the feathers, with further mass 
spectrometry analysis confirming that these pigments were 
carotenoids. This conclusion was confirmed by searching 
molecular weights in databases (METACYC, LIPID, KEGG), 
with a focus on 40 specific carbon compounds. According to the 
common carotenoids documented in birds (Hill & McGraw, 
2006), we suggest that the identified chemicals might be lutein, 
zeaxanthin, and echinenone. Therefore, it follows that at least 
two common dietary carotenoids (lutein and zeathanxin) exist in 
golden pheasant plumage. 

As the plasma samples achieved the same results in 
carotenoid identification, we inferred that this kind of dietary- 
derived pigment was absorbed and transported into feather 
follicles. Consequently, we determined the expression of 
several carotenoid-binding genes in golden pheasant feather 
follicles and found some patterns were significantly different 
between carotenoid-containing samples and mantle samples, in 
which carotenoids were absent in barbs. 

Genes related to carotenoid deposition did not show 
significant differential expression through RNA-seq analysis. 
Consequently, we performed real-time PCR and found that four 
genes (SfAR4, GSTA2, Scarb1, and APOD) were highly 
expressed in red breast and orange nape feathers, suggesting 
their possible function in carotenoid deposition in these tissues 
of golden pheasant. The differences in expression indicate that 
these genes might play different roles in carotenoid coloration in 


different feathers. In addition, all candidate genes were 
expressed in all tissues, even iridescent mantle feather follicles, 
while, the expression levels in yellow rump feathers showed no 
differences compared with iridescent mantle feathers. To clarify 
this issue, we tested a mantle feather using Raman 
spectroscopy and found carotenoids in the rachis, but not in the 
barbs. This led to the hypothesis that the transported 
carotenoids might combine with some other factors in the 
mantle feather rachis, such as keratins (McGraw et al., 2003), 
and thus the expression of carotenoid-binding genes might not 
be the only reason for carotenoid deposition in these regions. 
Compared with the mantle feathers, GSTA2 and APOD were 
highly expressed in red breast and orange nape feathers, but 
not in yellow rump feathers in which carotenoids were detected. 
Therefore, a further comparison with white control feathers was 
carried out. Results indicated that the two genes might be 
important to carotenoid deposition in golden pheasant plumage. 

The genes involved in carotenoid deposition in feathers are 
poorly known, yet are of substantial interest. The golden 
pheasant is an excellent species to address these questions, 
being one of the few Galliformes to express carotenoids in 
feather follicles. Six homologous candidate genes documented 
and investigated in Quelea quelea (Walsh et al., 2012) were 
studied in the golden pheasant in the present study. We found 
at least two candidate genes associated with carotenoid 
coloration, with differences in expressions raising the possibility 
that these genes might play different roles in golden pheasant 
plumage. Whether functional features of potential genes are 
species specific in golden pheasant or universal across avian 
taxa will spur the exploration of genetic and genomic 
investigations in future work. 
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